The Impact of Non-zero 6> 13 on 
Neutrino Mass and Leptogenesis in a SUSY SO (10) Model 
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The recent measurement of the reactor angle as sin 2 2#i3 = 0.092 ± 0.016(stat) ±0.005(syst) come 
from the Daya Bay collaboration. Evidence of nonzero #13 was also there at T2K, MINOS and 
Double Chooz experiments. We study the implication of these recent data on neutrino mass matrix 
and consequently on leptogenesis in a supersymmetric SO (10) model. To explain the smallness of 
neutrino mass, in general, we require a heavy Majorana neutrino which is a natural candidate in 
SO(10) model. In minimal SO(10) model, the symmetry breaking scale or the right-handed neutrino 
mass scale is close to the GUT scale. It is not only beyond the reach of any present or future collider 
search but the lepton asymmetry generated from its decay is in conflict with the gravitino constraint 
as well as unable to fit the neutrino data. We show that addition of an extra fermion singlet can 
accommodate the observed recent neutrino data in a supersymmetric SO (10) model. This model 
can generate the desired lepton asymmetry and provide TeV scale doubly-charged Higgs scalars to 
be detected at LHC. 
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The standard model (SM) of particle physics is based 
on the gauge group SU(2>) C x SU{2) L x U(1) Y , where 
C, L and Y respectively stand for colour, left-handed 
and hypercharge quantum numbers. In the SM, due 
to conservation of lepton family number and absence of 
any right-handed (RH) counterpart, left-handed neutri- 
nos are massless. However, it is well established that 
neutrino flavour oscillates which require neutrinos to be 
massive. One can generate tiny neutrino mass either by 
various type of seesaw mechanisms [1] or via loop cor- 
rections [2] going beyond the standard model. Neutrino 
oscillations can be parametrized using two mass squared 
differences, three mixing angles and one CP-phase of 
Pontecorvo-Maki-Nakagawa-Sakata matrix [3]. The third 
mixing angle, namely the reactor angle, was little known 
till date. 

In neutrino physics, a breakthrough measurement of the 
third mixing angle 813 come from Daya Bay experiment 
[4], and is confirmed by RENO experiment [5]. A more 
than 5cr measurement given by Daya Bay as sin 2 2^13 = 
0.092 ± 0.016(stat) ± 0.005(syst) and the corresponding 
value from the RENO experiment is sin 2 2#i3 = 0.113 ± 
0.013(stat) ± 0.019(syst). Their central values are very 
close to what predicted in ref. [6, 7] using the combined 
data set of T2K [8] and MINOS [9], earlier in the mid of 
last year, with more than 3a evidence. There was also a 
similar result from Double Chooz experiment [10] as well. 
All these new data, thus, make the neutrino mass matrix 
much more constrained and their consequences to other 
areas of physics. This new measurement led to prediction 
and implication of the # 13 angle in different ways [11]. 

The main ingredient of see-saw mechanism is heavy Ma- 
jorana neutrino, which is a natural candidate in the left- 
right symmetric SO(10) grand unified theories (GUTs) 
[12] . Embedding supersymmetry (SUSY) in such a model 
have some good features like protecting the Higgs mass 
from radiative correction that appear due to the huge 



difference between the weak and unification scales and 
help to have a good unification of gauge couplings at 
the GUT scale. We, thus, work here in a scenario of 
left-right symmetric (LR) supersymmetric SO(10) GUT 
model. In order to keep the gauge coupling unification 
intact at the GUT scale, any intermediate scale, here the 
LR-symmctry breaking scale, has to be very close to the 
GUT-scale. However, such a heavy symmetry breaking 
scale or an equally heavy Majorana neutrino is beyond 
the reach of any present or future collider analysis as 
well as is unable to produce the low energy light neu- 
trino data. 

Observed baryon asymmetry of the universe is another 
interesting problem. A popular explanation is to gener- 
ate baryon asymmetry via sphaleron process from lepton 
asymmetry [13, 14]. The later generally can be produced 
through the C and CP-violating out-of-equilibrium de- 
cay of heavy Majorana neutrinos, which is a member of 
SO(10)-GUT model and also responsible to explain the 
tiny neutrino mass. 

In the standard thermal leptogenesis, with heavy hi- 
erarchical right-handed neutrino spectrum, the CP- 
asymmetry and the mass of the lightest right-handed 
Majorana neutrino are correlated. In order to have the 
correct order of light neutrino mass-squared differences, 
there is a lower bound on the mass of the right-handed 
neutrino, Mjv^ 10 9 GeV [15], which implies a reheating 
temperature JllO 9 GeV. This will lead to an excessive 
gravitino production and conflicts with the observed data 
as discussed below. 

Gravitino, being the lightest and stable, is a suitable 
dark matter candiate in a R-parity conserving SUSY. In 
the post-inflation era, these gravitino are produced in a 
thermal bath due to annihilation or scattering processes 
of different standard particles. The relic abundance of 
gravitino is proportional to the reheating temperature of 
the thermal bath. One can have the right order of relic 
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dark matter abundance only if the reheating temperature 
is bounded to below 10 7 GcV [16, 17]. 

In this article we work in left-right symmetric SUSY 
SO(10) GUT model, rich with an extra SO(10) singlet 
lepton per generation [18-20] . This extra singlet lepton is 
a natural member in Eq and many other models. The is- 
sue is earlier addressed in different context [21, 22]. How- 
ever, here, we are to accommodated the recent neutrino 
data in this model. In addition, we discuss the impact 
of these new data on other related phenomenology. Our 
analysis, in a single model, is able to explain various is- 
sues like the light neutrino masses and their mixing as 
measured in recent experiments, have an exact unifica- 
tion of different gauge couplings at the GUT scale, have 
a low intermediate scale or a lighter right-handed Ma- 
jorana neutrino as well as to generate right amount of 
lepton asymmetry to explain the observed baryon asym- 
metry of the universe without being in conflict with the 
gravitino constraint. 

The spontaneous symmetry breaking prescription of 
SO(10) group in our model is as follows - At the GUT 
scale SO(10) is spontaneously broken with a simulta- 
neous vacuum expectation value (vcv) to the 3> 54 of 
along the direction singlet under the Pati-Salam group 
(Gps) SU(2) l x SU(2) R x SU(4) C [23] and to the sin- 
glet direction under the left-right gauge group (Glr) 
SU(2) L x SU(2) R x U(1) (B -L) x SU(3) C in the G PS 
multiplet (1,1,15) contained in a of SO(10). At 

this stage D-parity remains intact and both the gauge 
couplings of SU(2)l and SU(2)r. are equal, g^ = gn 
[24]. At the next step a vev to the D-Parity odd singlet, 
also contained in 4^2)° of SO(10), breaks the D-parity. 
To break the LR-symmetry at the next step we assign vev 
to the RH doublets xr © x R C 16h © 16h, however the 
subtlety of this breaking will be discussed later in the 
sections while at the last step electroweak symmetry is 
broken by a lOn-plet. 

The effective Lagrangian at the intermediate, LR- 
symmetry breaking, scale is 

Cy = r^ 6 ^«$10 + ^16^16^6 

+F^fT 1 X 1 R + fiT^T 1 + E.c. (1) 

The interacting superpotential to the scalar fields at the 
intermediate scale is given by 



w= M Ar A r ™A r ™ + M XrX r ° Xr ° 

+AAre^ + A'Af« (2) 

where ipL/R are left- /right-handed lepton doublets, while 
the superscript 16 stands to represent that they belong 
to the 16-plet of SO(10) representation and so on, and 
T, not to be confused with the superscript T for the 



transposed field, the fermion singlet field, one for each 
generation. The introduction of the scalar field A C 126, 
we can justify from the scalar interaction terms as follows. 

We have assigned a vev to the right-handed doublet com- 
ponent only {Xr 6 ) = (Xr 6 ) = v x- However, the vanish- 
ing F-term conditions give us, 

(A- 60 ) = vr = —A (3) 
M Ar 

We, thus, have a large induced vev to the neutral com- 
ponent of the triplet scalar A R or A R , once the neutral 
doublet component Xr gets a vev. For example, with a 
lighter RH-triplct mass M A ~ 100 GeV- 1 TeV, it is pos- 
sible to have v R ~ 10 9 -10 12 GcV for v x = 10 6 -10 7 GeV, 
assuming A ~ 0(1). Since, vr 3> v x , the spontaneous 
symmetry breaking of the group SU(2) R x U(1)b-l — > 
U(l)y takes place at a higher scale generating large RH 
Majorana neutrino masses Mn 3> Mx- This will lead to 
a small N% — Tj mixings, which is a crucial point needed 
to establish the out-of equilibrium conditions for lepto- 
genesis. 

In this model the neutral fermions per generation are a 
left-handed neutrino is, a right-handed neutrino N, both 
of which are member of 16-plet of SO(10), and a sterile 
neutrino, T. From the Yukawa interaction, Eq.(l), we 
see in the (u, N, T) basis the 3x3 mass matrix is given 
by 

/ m D \ (u c \ 
M v = (v N c T ) L [ml M N M x V . (4) 
V Ml M J \ tJ l 

Here the N — T mixing matrix arises through the vev of 
the RH-doublet field with 

M x = Fv x , where v x = ( X ° R ), (5) 

and the RH-Majorana neutrino mass is generated by the 
induced vev of the RH-triplet with 

M N = fv R , (6) 

where, vr is given in eq.(3) as described above. The 
vev of the weak bi-doublct $(2, 2, 0, 1) C 10 H of SO(10) 
yields the Dirac mass matrix for neutrinos, 

m D = y($°>. (7) 

In our model different mass scales hierarchy is Mjy 3> 
Mx ^> M ^ m D- With this hierarchical mass spectrum 
[20] , the eigenvalues of the mass matrix in eq. (4) is given 
by 

m v m D [M^ 1 (J* {Ml)- 1 ] m T D , (8) 
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Here, we see that the light neutrino masses satisfy a dou- 
ble see-saw structure. It may be noted that the mass ma- 
trix structure in eq. (4) ensures that the type-I see-saw 
contribution is absent and Mx remains unconstrained 
by the light neutrino masses. This freedom in Mx - a 
hallmark of the model - is vital to ensure adequate lep- 
togenesis. 

In order to satisfy both the neutrino data as well as to 
generate the required amount of lepton asymmetry in this 
SO(10) model, we note that the mass matrix /i can be 
obtained using eq.(8) once we know the mass matrices 
m„, mrj and mx- Our strategy is described as follows: 

To construct the Dirac mass matrix mn, here we work in 
a basis in which the down-quark and charged lepton mass 
matrices are diagonal. The entire mixings in the quark 
and lepton sectors, thus, can be ascribed to the mass 
matrices of the up-type quarks and the neutrinos, respec- 
tively. On the otherhand, in SO(10) model, the quark- 
lepton symmetry [23] relates the neutrino Dirac mass ma- 
trix mo to its counterpart in the up-quark sector. We, 
therefore, obtain mo using the quark masses and the 
Cabibbo-Kobayashi-Maskawa mixing angles, upto 0(1) 
effects due to RG evolution. Using the data enlisted in 
the Particle Data Group [25] for the CKM matrix ele- 
ments, its Dirac phase, and the running masses of the 
three up-type quarks, namely, m u = 2.5 MeV, m c = 1.29 
GeV, m t = 172.9 GeV, we have, 

m D ^ M v = Vq KM diag(m u ,m c ,m t )VcKM, (10) 

where we have used the CKM phase 6c km = 1.2 radian, 
and the quark mixing angles sm8^ 2 = 0.2253, sinfl^ = 
0.041, and sin#^ 3 = 0.0041. The Dirac neutrino mass 
matrix is fixed by the underlying quark-lepton symmetry 
of S'O(IO). Neglecting small RG corrections, it is taken 
to be approximately equal to the up-quark mass matrix. 

Next, we see that the matrix Mx is determined through 
eq.(5). However, the 3x3 coupling matrix F is 
completely arbitrary. To minimize the number of in- 
dependent parameters, we take the matrix F to be 
real and diagonal. Here, we choose, for example, 
M x = diag(0.15, 0.5, 0.8) x v x with v x = 10 6 GeV. Once 
a vev to the RH doublet xjt is chosen, we can have a 
definite induced vev to the RH triplet Aj^ 26 using eq.(3). 
Not to mention, due to same reason like Mx, we choose 
/ <~ diag(0.1, 0.5, 0.9) to obtain, via eq.(6), the Mx mass 
matrix of O(10 10 ) GeV. This will, for a general M x , lead 
to a Ni — Tj mixing given by 



sin6*i; 
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M Ni 



(11) 





Sm 2 
/10" 5 eV 2 


Am 2 
/10" 3 eV 2 


sin 2 9i2 


sin 2 023 


sm 2 6li3 
/lO" 2 


bfv 


7.58 


2.35 


0.312 


0.42 


2.5 


Iff 


7.32 - 7.80 


2.26 - 2.47 


0.296 - 0.329 


0.39 - 0.50 


1.8 - 3.2 


2a 


7.16 - 7.99 


2.17 - 2.57 


0.280 - 0.347 


0.36 - 0.60 


1.2 - 4.1 


3a 


6.99 - 8.18 


2.06 - 2.67 


0.265 - 0.364 


0.34 - 0.64 


0.5 - 5.0 



TABLE I: Ranges for mixing parameters obtained in Ref. [6] 



In neutrino physics, breakthrough measurement of the 
third mixing angle #13 come from different experiments. 
Its evidence come in the mid last year from T2K and 
MINOS. Their combined data predicted non-zero #13 [6, 
7] with more than 3(7 evidence. Recently, Daya Bay as 
well as RENO experiment come up with a more than 5cr 
measurement. However, their central value is very close 
to what predicted in [6, 7]. To construct the neutrino 
mass matrix we use the combined neutrino mixing data, 
including the recent T2K and MINOS results, given in 
Table-I from Ref. [6] with the new reactor flux estimate. 
In addition, here, we assume all other CP-phases in the 
lepton sector, except the one through the CKM matrix 
in the quark sector, to be zero. Assuming the lightest 
neutrino mass eigenvalue (mi for normal hierarchy and 
TO3 for the inverted hierarchy) to be zero, we obtained two 
other mass eigenvalues using different values of 5m 2 (= 



m{) and Am (= 



(m| + ml)/ 2) from Table- 



I. We construct m u from these mass eigenvalues via the 
PMNS matrix, Upmns as 



m v = Up MNS diag(mi,77i2,m 3 )[/pA/jvs 



(12) 



With the knowledge of m„, mry and Mx we then use 
the inverse see-saw mass formula, from eq.(8), to obtain 
elements of the matrix \i for both normal and inverted 
hierarchical light neutrino masses. The Mx and Mx 
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FIG. 1: Normal Hierarchy: T\ mass corresponding to the best 
fit value (cyan-cross) of neutrino data and its variation for 
la, 2a and 3<r are shown by the red (central dark) , yellow 
(whitish) and brown (outer dark) areas. 
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3 3<J range 




FIG. 2: Inverted Hierarchy: Ti mass corresponding to the 
best fit value (cyan-cross) of neutrino data and its variation 
for la, la and 3a are shown by the red (central dark) , yellow 
(whitish) and brown (outer dark) areas. 



matrices are used in eq. (9) to compute the mass eigen- 
values of the singlet fermions and their mixings. Thus 
the two input matrices M^y and Mj (chosen diagonal) - 
eq. (4) - completely determine the singlet neutrino, T i: 
masses and mixings consistent with the recent data on 
the light neutrino mass spectrum and their mixing. Out 
of three eigenvalues only one, denoting it as Ti, is above 
the threshold energy to decay into l(f>. We have shown 
the variation of the allowed T\ mass in Fig. 1 for a nor- 
mal hierarchical light neutrino mass and in Fig. 2 for the 
case of inverted hierarchy. For normal hierarchy, best fit 
values given in Table -I corresponds to a T\ mass equals 
to 0.6124 x 10 4 GeV and is denoted by the cross (cyan) 
in Fig. 1. Here, we see that T\ mass varies between 
(0.4709 - 0.8526) x 10 4 GeV (brown-central dark) for la, 
0.2874 x 10 4 - 0.1152 x 10 5 GeV (yellow - whitish) for 2a 
and 0.1388 x 10 4 - 0.1533 x 10 5 GeV (red - outer dark) for 
3cr allowed experimental neutrino data. A similar analy- 
sis is shown for the inverted hierarchy is given in Fig. 2. 
With the same Mx and Mn, for inverted mass hierarchy 
light neutrino case, the best fit values of Table -I corre- 
sponds to a Ti mass equals to 0.1585 x 10 6 GeV and is 
denoted by the cross (cyan) in Fig. 2. In this case, Ti 
mass varies between (0.1556 — 0.1624) x 10 6 GeV (brown- 
central dark) for la, (0.1504-0.1662) x 10 6 GeV (yellow 
-whitish) for 2a and (0.1427 - 0.1713) x 10 6 GeV (red - 
outer dark) for 3cr allowed experimental neutrino data. 

We now discuss to check if or not the predicted mass 
spectrum for Ti, obtained using the neutrino data, is 
able to generate require amount of lepton asymmetry. 
The mass scale Ti is well-below the condition on reheat- 
ing temperature come from the gravitino overproduction. 
This model, thus, is in good agreement with the current 
limit on the dark matter relic abundance. 

The singlet fermions decay through their mixing, con- 
trolled by the ratio Mx/Mm, with the Ni. The latter, 



FIG. 3: The tree and one-loop contributions to the decay of 
Ti that generates the lepton asymmetry. 



which have masses O(10 10 ) GeV and are off-shell, decay 
to a final l<f> state, where I is a lepton doublet and <j) the 
standard Higgs boson. This two-step process - for which 
a typical tree diagram is depicted in Fig. 3 - results in 
a lepton asymmetry of the correct order. Because of the 
large value of M n >> Mx , a small Tj — JVf mixing is nat- 
urally permitted which in turn guarantees out-of- equilib- 
rium condition to be realised near temperature T ~ Mt- 

Below we discuss various solutions of the Boltzmann 
equations. These determine the number densities in a 
co- moving volume Yr = tit j rig and Yt, = tit, Ins, where 
riL and ns are respectively the number densities of lep- 
tons and the entropy. We can read the equations as - 



dY T 
dz 
dY L 

dz 



(Y T 



er TTf.\ (Xt 



zH(z) 



1 D 

zH(z) 



zH(z) 

rf 



zH(z) 



(13) 



where T^, and T^y represent the decay, scattering, 
and wash out rates, respectively, that take part in es- 
tablishing a net lepton asymmetry. We refrain from pre- 
senting their detailed expressions here and due to neg- 
ligible contribution from supersymmetric processes [26], 
we have not included them. The Hubble expansion rate 
H(z), where z = Mt/T, and the CP-violation parameter 





10 100 



FIG. 4: The decay and inverse-decay rate of T are compared 
with the Hubble expansion rate, H, as a function of z, for 
the best fit values of neutrino data only are shown for the 
Normal Hierarchy (Left) and Inverted Hierarchy (Right). 
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FIG. 5: Normal Hierarchy: The comoving density of T 
- and the leptonic asymmetry - Yl 



Y t 



as a function of z is 
shown. Corresponding to the best fit neutrino data, Yl, Yt 
and Y^(® are shown in the central (cyan) lines and the dot- 
dashed (red), dashed (maroon) and solid (blue) boundary lines 
corresponds to la, la and 3a ranges. 



are given by 

H(M T ) 1/2 M| 
H(z) = - — , H(AIt) = 1.67.9*' 



M„ 



r(T -> i<f>) - r(T -> i(j>* 

T(T -> ld>) + r(T ->• l(j)* 



(14) 



In Fig. -4, we have shown the variation of the decay and 
inverse decay rate of T\ for both normal (left) and in- 
verted (right) hierarchical light neutrino cases. Corre- 
sponding to the best fit value mass spectrum, the fig- 
ure clearly shows how the out of equilibrium condition 
are satisfied to generate the lepton asymmetry. We have 
shown the lepton asymmetry production results in Fig. -5 
and in Fig. -6. We assume that in the very initial stages 
the number densities Tj, i = 1,2,3, are zero. T\ decay 
through the channel l<f> to produce the lepton asymmetry. 
One important point to note here is that in this process 
of lcptogcnesis, reheating temperature is consistent with 
the gravitino constraint. In the figure, corresponding to 
the best fit values of neutrino data, Yl, Yt x and V^? 
are plotted as the central (cyan) lines. The effect of la, 
2(7 and 3cr ranges of neutrino data on the evolution of 
Yr and Yr x are shown respectively with the dot-dashed 
(red), dashed (maroon) and solid (blue) boundary lines. 
It is seen from the figure that although it is perturbed at 
a lower value of z but as the universe expands Yr achieves 
the right order (~ 10~ 10 ) starting off from a vanishing 
initial value while that for Yr t are well separated. How- 
ever, for the inverted hierarchy case, effect for la, 2a 
and 3a are overlapping due to a relatively monochro- 



FIG. 6: Inverted Hierarchy: Same as in Fig. 5 



matic and large T\ mass. Now, we just comment on how 
to achieve a unification of gauge couplings in this SO(10) 
model. As mentioned earlier, Higgs multiplets 210 and 
54 arc utilised to break the symmetry at Mjj. Within 
the 210 there are two components which develop vevs; 
one breaks SO(10) to C/3221 while the other is respon- 
sible for D-parity breaking. The SU(2) R x U(1) B -l 
symmetry is broken by the induced vev ~ 10 11 GeV, 
which is also responsible for the masses of the iVj, of RH 
triplets in 126 © 126. The last step of breaking car- 
ried out by the weak bi-doublet in 10. With the anal- 
ysis of the gauge couplings RG evolution we determine 
the intermediate mass scales. An intermediate scale at 
Mr ~ 10 9-11 GeV can be obtained through the intro- 
duction of effective dim. 5 operators scaled by the Planck 
mass, Mpi [27]. It is interesting to note that both 210 
and 54 are necessary for viable SUSY SO(10) breaking 
pattern and consequently the resulting two dim. 5 opera- 
tors appear to alleviate the problem of leptogenesis under 
gravitino constraint: 



£nro 



2M, 



-Tr [F m „(t?i$ 2 io + m^i)F^] . (15) 



The reason behind this is that above interaction lead to 
finite corrections to the gauge couplings at the GUT- 
scale so that the gauge couplings of left-right gauge group 
emerge from one effective GUT-gaugc coupling. The up- 
shot of this is that with these additional contributions it 
is possible to lower Mr to as low as 10 9 — > 10 11 GeV 
as required in this model. The grand unification scale is 
as large as M\j ~ 10 17-18 GeV and the model predicts 
a stable proton for all practical purposes. Another way 
to achieve this gauge coupling unification is to introduce 
some additional scalar multiplet at the intermediate scale 
[22]. 
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Finally, wc comment on the experimental prospect of 
doubly charged scalar of this model at LHC or ILC [22]. 
Using the D-parity mechanism in this model we make 
the RH-triplets in 126 H ® 126 H carrying B — L = ±2. 
Their masses are from 100 GeV to a few TeV. This light 
triplet scalar comes out as a necessary condition to en- 
hance the induced vev, vr or to break the LR gauge sym- 
metry at high scale. Consequently, we have heavy RH 
Majorana neutrinos as well as the corresponding gauge 
bosons. This forbids to decay into right-handed 

gauge bosons. However, after being produced via Drcll- 
Yann process at LHC, these doubly-charged scalars will 
decay to fermions to be detected at LHC. 

In conclusion, in view of the recent neutrino data we have 
presented a left-right symmetric SUSY SO(10) model. 
This model is capable to solve a multi-dimensional prob- 
lems. This model has the following features 

• By virtue of its construction, it is consistent with 
the most recent neutrino masses and mixing an- 
gles obtained at MINOS, T2K, Daya Bay, RENO 
experiments. 

• We have discussed the variation of the T\ mass due 
to lcr, 2cr, 3cr variation for both the normal and 
inverted hierarchy light neutrino cases. 

• It generates a correct lepton asymmetry via the de- 
cays of SO(10) singlet neutrino with a mass scale 
to be consistent with the gravitino constraint. 

• It can also have a good unification of different gauge 
couplings at the GUT scale. 

• The model is also a source of light doubly-charged 
scalar. It's mass range is within the reach of the 
LHC. 
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